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Abstract 

Self-assembled quantum dots (QDs) are highly strained heterostructures, the lattice strain signifi¬ 
cantly modifies the electronic and optical properties of these devices. A universal behavior is observed in 
atomistic strain simulations (in terms of both strain magnitude and profile) of QDs with different shapes 
and materials. In this paper, this universal behavior is investigated by atomistic as well as analytic con¬ 
tinuum models. Atomistic strain simulations are very accurate but computationally expensive. On the 
other hand, analytic continuum solutions are based onassumptions that significantly reduce the accuracy 
of the strain calculations, but are very fast. Both techniques indicate that the strain depends on the 
aspect ratio (AR) of the QDs, and not on the individual dimensions. Thus simple closed form equations 
are introduced which directly provide the atomistic strain values inside the QD as a function of the AR 
and the material parameters. Moreover, the conduction and valence band edges Ecjv and their effective 
masses rn*Qiy of the QDs are dictated by the strain and AR consequently. The universal dependence of 
atomistic strain on the AR is useful in many ways; Not only does it reduce the computational cost of 
atomistic simulations significantly, but it also provides information about the optical transitions of QDs 
given the knowledge of Eq/v and tuq/y from AR. Finally, these expressions are used to calculate optical 
transition wavelengths in InAs/GaAs QDs and the results agree well with experimental measurements 
and atomistic simulations. 

Keywords — Self-assembled quantum dots, Stranski-Krastanov, Atomistic strain. Analytical continuum 
strain. Continuum elasticity. Optical transition. 


1 Introduction 


Self-assembled quantum dots (QDs) have improved the performance of many optoelectronic devices, such as 
quantum dot infrared photodetectors (QDIPs) [^[^, intermediate band solar cells (IBSCs) [^|^, quantum 
dot optical amplifiers [^, and quantum dot lasers [^[^; QDIPs have lower dark current than conventional 
photodetectors [2l^ and are sensitive to normally incident light unlike their counterparts that are made 


from quantum wells 2 9 10 . For IBSCs, one of the most successful methods in pushing solar cell efficiency 


beyond the Shockley-Queisser limit is to add one or more intermediate bands inside the gap, which can be 
realized by using quantum dots and quantum dot-in-a-well devices [^. Quantum dot optical amplifiers offer 
ultra wide band polarization insensitive high-power amplification [^. Quantum dot lasers provide the most 
reliable and temperature-insensitive operation around the spectral range of 1.15 to 1.25 ym [^. 
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Self-assembled quantum dots are highly strained heterostructures. The atomistic strain in such structures 


is usually on the order of 10%. The electronic band structure is affected significantly by the lattice strain 11 


For proper band structure calculations, the strain needs to be calculated and included accurately 12 . Such 


high strain values are beyond the domain of validity of the continuum elasticity theory [13| and a rigorous 
atomistic treatment of strain is needed (such as Keating 14 or anharmonic models 11 1^). 


Some previous works have calculated strain in QDs using continuum elasticity theory [W -18 . These 


continuum elasticity calculations become inaccurate as the strain exceeds the elastic limit and goes beyond the 
domain of validity of continuum elasticity, which is the case in these self-assembled QDs. TTp3 . In addition, 
the mathematical modeling problem is complicated and many assumptions are needed to obtain an analytic 
solution. These assumptions further reduce the accuracy of the analytic solution. Here, these assumptions 
are discussed and avoided as much as possible for a better accuracy. Nevertheless for a precise treatment 
of strain, atomistic simulations need to be used. The main issue with atomistic simulations is that they are 
computationally very expensive. In this work, both approaches are investigated and compared against each 
other. The analytic solution provides qualitatively similar results to the atomistic simulations. However, it 
quantitatively underestimates the strain magnitudes. On the other hand, the atomistic simulations of QDs 
show a universal behavior in the strain results, independent of QD dimensions and materials. Accordingly, 
compact expressions which capture this universal behavior are put forward as replacements of the analytic 
solution. These closed form expressions provide the atomistic strain values without the numerical burden of 
atomistic simulations and the inaccuracy of analytic models. 

In section H, the analytic continuum solution is discussed. Later on, the atomistic results are shown in 
section HI. Finally, the universal behavior is introduced and its applications are discussed in section IV. 


2 Analytical Continuum Solution 


In this section, analytic expressions for the strain in cuboid QDs are introduced. These expressions are 
obtained by solving the continuum elasticity equations [T^ which are based on the stress field Green’s function 


17 19 20 used to predict stain in quantum wells, wires, and dots. Unlike previous works 16 17 19 20 


hydrostatic strain is not assumed to be a constant here. The assumptions and details of the derivations are 
provided in the appendix. The strain in the middle of a cuboid quantum dot with base length b and height 
h is found to be: 


£h 


l-2v 

E 


127rA -I- 16Atan ^ 




( 1 ) 



where Eh is the hydrostatic strain, eh = Exx + Syy + £221 and £_b is the biaxial strain, Eb = Exx + £yy — ‘^Szz- 
A = . ' , where En is the lattice misfit, En = ^ gj.g pgisson ratio and 

Young’s modulus of the quantum dot material. The choice of describing the strain in terms of Eh and Eb 
components instead of Cartesian components £yy and Ezz is due to the fact that in zinc-blend materials 
the deformation in band edges depends on the Eh and Eb directly. 

The analytic continuum strain solution indicates that the strain depends on the base to height aspect 
ratio (b/h). Figureshows Eh and Eb as a function of b/h for a cuboid QD obtained from equations Q 
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and ([^. The magnitude of Eh decreases slightly with increasing b/h, while the magnitude of es increases 
with increasing b/h. 



Figure 1: The analytic continuum solution of hydrostatic and biaxial strain as functions of b/h. Strain 
depends on b/h, but not on the individual dimensions. 


Notice that the solution in is based on a stress field Green’s function that has been introduced for 
in which at least one strain component is assumed to be known. This assumption is 


20 


quantum wires 

suitable for quantum wires and quantum wells but not for quantum dots. The main reason behind the 
deviation from the atomistic simulations of the QDs is this assumption. It is shown in the next section that 
the analytic solution converges to the atomistic results at large b/h values since the QD becomes more like 
a quantum well in which this assumption is more justified. 


3 Atomistic Simulation. 


For an accurate calculation of the strain, a rigorous atomistic strain model is needed to calculate the relaxed 


atom positions. The Keating model 14 provides the total elastic strain energy of the system as a function 
of atom positions. 
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(3) 


where the coefficient a corresponds to the force constant for the bond length distortion, and /3 corresponds 
to the bond angle distortion as shown in Figure a and /3 are material constants, r^n is a vector from 
atom m to atom n for the strained crystal, while dmn is the same vector for the bulk unstrained crystal. 
These atomistic interactions are assumed to be between the nearest neighbors only. 


Keating, in his seminal work 14 , showed that there is a straightforward connection between the atomistic 


force constants a and (3 and the macroscopic elastic constants (e.g. cn and C 12 ) of the material. The Keating 
model is known to be a suitable model for calculation of the atomistic strain in self-assembled QDs Tl|[^ - 23 


and phonon dispersion in nanowires and in bulk 24 . 
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Figure 2: Bond lengths and bond angle for three neighbor atoms m, n, and k. 


Cuboid and dome shaped QDs with different dimensions have been simulated using the Keating model 
to study the behavior of strain with different QD dimensions, shapes, and materials. The dimensions of the 
whole simulated domain are 60 nm x 60 nm x 60 nm and the quantum dot is located in the middle of the 
structure. The strain simulation contains about 10 million atoms with the atomistic grid shown in Figure 
Such large systems are computationally expensive for strain simulations and require a highly scalable 


parallel code. All the atomistic strain simulations in this work have been done by NEM05 25 26 



Figure 3: The atomistic grid of a cuboid InAs/GaAs QD with the base length of 24 nm and height of 4 nm. 

Figure shows atomistic Eh and Eb along a line in the middle of the device in [001] direction. The 
positions of the atoms on the line are divided by the dot height (h) for comparison purposes. QDs with 
the different dimensions but the same aspect ratio have similar strain distribution over the normalized 
coordinates. This behavior is observed in both dome and cuboid QDs. 
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InAs/GaAs Cuboid QD InAs/GaAs Dome QD 



z/h z/h 

Figure 4: Atomistic Eh (blue) and Eb (red) strain profiles along a line in the middle of the QD with the heights 
of 2nm (lines) and 4nm (symbols) in [001] direction. The solid lines and symbols overlay each others. QDs 
with different dimensions but same aspect ratios exhibit same strain profiles over the normalized coordinate 
(z/h) for different aspect ratios (a)-(f). The left three figures are strain distributions of cuboid QDs with the 
base to height ratios (b/h) of 1 (b), 2 (b), and 5 (c), while the right three figures are for the dome shaped 
QDs with the base diameter to height ratios (D/h) of 1 (d), 2 (e), and 5 (f). The main point is that, QDs 
with certain shape and aspect ratio have a unique strain distribution. 

Figure depicts the atomistic Eh and Eb versus the analytic solution in the middle of a cuboid quantum 
dot. The atomistic results show similar behavior to the analytic solution qualitatively but with significantly 
higher magnitudes. Similar trends exist in the other III-V material systems and shapes, as shown in Figs, 
and respectively. The universal behavior that the atomistic strain in different materials depends on the 
aspect ratio of the QDs independent of actual physical dimensions can be captured with compact equations 
and be used to approximate strain in any QD, as discussed in the next section. 
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Figure 5: a) The atomistic and analytic Eh and Eb at different base to hight ratio (b/h) of cuboid InAs/GaAs 
QD. The analytic solution underestimates the strain magnitudes. The atomistic Eh and Eb with different 
materials in cuboid (b) and dome shaped (c) QDs as a function of aspect ratios. The strain depends on the 
aspect ratio (b/h or D/h) and not on individual dimensions. 


4 Universal Behavior and its Applications 

As described in the previous section, atomistic strain depends only on the material properties and the aspect 
ratio of the QD. This behavior can be captured with compact equations as follows; First a function of aspect 
ratio is defined. This function is called universal behavior function U{AR), where the aspect ratio AR is 
dehned as AR = ^ in cuboid QDs (the definition of AR in dome shaped QDs is discussed later). The strain 
in the middle of the QD can be written in terms of this function: 


e{AR) = e{AR ^ oo) + U{AR) [e{AR = 1) - e{AR oo)], 


(4) 


^ g(Aj?) - e{AR oo) 

^ ’ e{AR = 1) - e{AR ^ oo) 


(5) 


The universal function I/(Ai?) is different for hydrostatic (i.e. Uh {AR)) and biaxial strains (i.e. Ub {AR)), 
but both of them satisfy the following conditions: 


U{AR) 


0, AR —>■ oo 
1, AR=1 


( 6 ) 


Assuming that the strain values are known at the limits of AR = 1 and AR —)■ oo, the functions Uh (AR) 
and Ub {AR) can be fitted numerically to the atomistic simulation results using equation ([^ as follows: 


Uh {AR) 


1 

0.75 (Ai?)”° '‘+ 0.25 (Ai?)^ "^’ 


Ub {AR) 


1 

(Ai?)^-^®' 


(7) 

( 8 ) 


To obtain a universal behavior function, the QDs with different materials and shapes need to fit into one 
equation which is the case as shown in Figs. [6}r and [§r. 
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Figure 6: a) The hydrostatic Uh{AR) and b) biaxial Ub{AR) universal behavior data calculated for different 


materials in cuboid and dome shaped quantum dots. The black lines are 



obtained from equations Q and 


The strain values at the limits of AR = 1 and AR —>■ oo need to be determined still. The strain values are 
known for high aspect ratios {AR —>■ oo) since in that regime the QD becomes a quantum well with accurate 
analytic strain solution: er {AR oo) = — 2£o and Eb {AR —>• oo) = —2eqj^^, where £o is the lattice 

misfit, £o = , and Vd is the Poisson ratio of the dot’s material. On the other hand, in the case 

of AR = 1, there is no analytic solution available for the strain except for eb in cuboid QDs which is zero 
due to the symmetry: Eb {AR = 1) = e^x + Syy ~ 2£z2 = 0, the rest of the values are approximated. In the 
cuboid QDs, the atomistic simulations reveal that er is close to {AR = 1) = —2.77£o , where 

Vav is the average Poisson ratio of the dot and the substrate materials. Whereas in the dome shaped QDs, 
the atomistic strain results are close to = l) = —2.7£o and Eg°"^^ ~ ^) ~ 0.2£o . 

Combining the strain values at the limits of AR = 1 and AR —>■ oo with equation provides the strain 
equations in the cuboid and dome QDs: 
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( 10 ) 

( 11 ) 

( 12 ) 


The exact atomistic strain values in QDs with different materials show deviations less than 10% from the 
proposed universal function (shown in Fig. [^. The black solid lines are obtained from the equations ([^ to 
(12). Defining AR in the dome shaped QDs as AR = 0.85^^ makes it possible to use the same universal 


function in cuboid and dome shaped QDs. 

Equations 0 to ( |12[ ) can be used to estimate the atomistic strain in cuboid and dome shaped quantum 
dots with good accuracy. These expressions can guide researchers studying quantum dots to understand the 
trends of the atomistic strain in self assembled QDs. In addition, the optical and electronic properties of QDs 
depends mainly on the strain. The effective mass and band edges of QDs can be determined from the strain 
values. Moreover, a prior knowledge of strain inside QDs helps to design and simulate QDs with required 
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optical absorption spectrum. The effect of strain on the band edges can be calculated using deformation 
potential theory 27 : 

^Eq — H : 

^EyHH = dySH + 2 ^-®’ 

^EyLH = dySH — 2^B, 

where AEy is the shift in the conduction band edge due to the strain, AEyHH and AEy^n are the shifts in 
the heavy and light hole band edges, respectively, dy, dy, and b are the deformation potential coefficients of 
the material. 

To show the impact of these expressions, the optical transition energies of an InAs/GaAs dome shaped 
quantum dot with height of 5 nm and diameter of 20 nm are calculated here using different techniques and 

The effective mass simulation was performed based on 


compared against experimental measurements 28 


an efficient scheme . Including the universal strain behavior in the effective mass model (EM) reduces the 
error in the predicted optical transition from 25% to 4% as shown in Table Despite the error not being as 
small as a full band tight binding (TB) simulation (1.6%), the computational cost is much less as indicated 
in Table The tight binding simulations were performed in NEM05 using sp^cf’s* basis [^. The values of 
the InAs effective masses under strain are taken from |^. Figure]^ shows the optical transitions calculated 
from the effective mass model including the universal strain and the full band tight binding simulation. The 
deviation of the results between these methods does not exceed 4.5%. 



EM no Strain 

EM with Strain 

TB Atomistic 

Experimental 

Optical transition 

0.7169 eV 

0.9145 

0.9377 

0.976 - 0.93 eV 

Error (vs 0.953 eV) 

25% 

4% 

1.6% 


Number of CPUs 

1 

1 

216 


Time/CPU (s) 

80 

81 

7856 



Table 1: Calculated optical transition energy using different techniques versus the experimental measurement. 
Using the universal behavior equations to include the effect of strain on the band edges and effective masses 
significantly improves the effective mass model results without any additional computational cost. 













Figure 7: Optical transitions in InAs/GaAs QDs as a function of the dot height for different AR calculated 
by the effective mass model (symbols) including the universal strain versus the full band tight binding 
calculations (solid lines). The deviation does not exceed 4.5%. 


5 Conclusion 

Atomistic simulations are known to be the most accurate tool for studying QDs, however they are compu¬ 
tationally expensive. On the other hand, the analytic continuum solutions are very fast, but they are based 
on the assumptions that significantly reduce the accuracy of the predicted strain in QDs. Using the fact 
that atomistic strain in QDs depends on the aspect ratio of the QD {universal behavior) a novel method to 
predict the strain is proposed. It is shown that this universal behavior of atomistic strain exists in QDs with 
different shapes (e.g. cuboid and dome shaped QDs) and materials (e.g. InAs/GaAs, InP/GaP, etc.). The 
proposed compact strain equations can be used to design the optical properties of the QDs or predict the 
dimensions of a fabricated self assembled QD from its optical transitions with good accuracy. 
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Appendix: Analytical solution for strain in quantum dots 

In this section, we have derived simple analytical expressions that describe strain in self assembled quantum 
dots based on continuum elasticity theory. The derivation here uses the analytical solution for strain inside 
the quantum dot provided in [16|. The reader is advised to read 16 first. Before we derive the general 


expressions, it is helpful to understand the behavior of the strain in an ideal quantum well first. 


Cuboid quantum dot 


In 16 , an analytical expression for the strain of a clipped pyramidal quantum dot has been obtained using 


the continuum elasticity and Green’s function technique. For notations and the meaning of the parameters 


used in the following derivation, please refer to the mentioned paper and its predecessors 16 17 19 20 


We will follow the same assumptions of the paper 16 except for assuming constant hydrostatic strain. 


Henceforth the assumptions of the proposed analytical solution are: 


• The solution is based on a stress field Green’s function that has been introduced for quantum wires 20 


for which at least one strain component should be known. This is suitable for quantum wires and 
quantum wells but not accurate for quantum dots. This assumption is the main reason behind deviation 
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from the atomistic result of the quantum dot as shown in section If. The matching of the analytical 
solution with the atomistic results at high values of base to hight ratio happens because the quantum 
dot becomes more like a quantum well. 

• Continuous, elastic, linear, and isotropic mediums. 

• Infinite substrate. 

• Elasticity constants of dot and substrate are the same. 

• The quantum dot is a cuboid with base dimensions a = b and height H. 

• The strain is calculated at a point in the middle of the quantum dot. 

The hydrostatic strain Eh is defined as. 


— ^xx T ^yy “b 
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(13) 


The biaxial strain es is defined as. 


E^ ^XX ~b Eyy 2e 


in a similar way to the hydrostatic strain. 


1 V , o \ 
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£b = <^b- 


(14) 


Following the expressions in 16 and using the assumptions mentioned before: The dot has a cuboid 
shape, 

/i —>00 , /—tO , E[ = hf^ 

since we are calculating the strain in the middle of the dot then, 

H 

Xl= X 2=0 , 0:3 = y. 


calculating the parameters in 16 under these assumptions with some algebraic simplifications yields, 

Xi = X 2 = X 3 = pa ,p = ±1, 

Y^=Y2 = Y3 = qb ,g = ±l, 

-H H 

^2-y- :^3-y, 

A = B = Hi = Q ,hHi = l, 


a^ + b^ + ^, 


82 = 83 = 3 = 
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Cii = 0 ^ + 62 + qbS 


Dll = paZn 


a = 


pa{S"^ + qhS) 




/? = 


pZn{b‘^ + qbS) 


P , ^+iH {b^+qbS) 
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(15) 

(16) 
(17) 


Since HI = 0 and hHl = 1 , the expression for an is simplified to only one term in the summation 

= -4^A + A y (-l)’^+igtan-i(-), (18) 

rv 


CTll 


p=±l 

q=±l 

n=2,3 


where A = substituting from equation (17) and simplifying the summation. 


an = —47rA + 4A ( tan 


_i ( H{b^ + bS) \ _i( H{b^ - bS) 


\2a{S‘^+bS) 


— tan 


\2a{S^-bS) 


(19) 


For a^x and ayy, put a = b which is the base length, and put H = h which is the height of the dot in 


equation (19) 


axx = o'yy = —47rA + AAtan 


-1 


2(!: 




— 4Atan 


.1 1 t-ym+1 


.hJ ' 4 h 


2{iY + \-lJ2{i) 


with some algebraic and trigonometric simplifications, one obtains 


, ( 20 ) 


axx = o'yy = —47rA + SAtan 


For azz, exchange a and H in equation (19) then put a = b and H = h, 


( 21 ) 


a^z = —47rA + AAtan 


v/ptM? 




— AAtan 


-1 




^ 1 +!«)'-! v'iTiot 

with some algebraic and trigonometric simplifications, one obtains 


• ( 22 ) 


13 



























47rA + 8 Acot ^ 


(23) 


Czz — — 



Substituting equations (211 and (23) into equations (|13| and (|14[), we get the final expressions for the 


hydrostatic and biaxial strain: 


£h 


1 - 2v 
E 


127rA + 16Atan ^ 




£b 


l + v 
E 


16A 





(24) 


(25) 


14 
















